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ABSTRACT 

The  local  theory  of  singular  points  is  extended  to  a  large  class  of 
linear,  second-order,  ordinary  differential  equations  which  can  be  physical 
Schroedinger  equations  or  govern  the  modulation  of  real  oscillators  or 
waves.  In  addition  to  Langer’s  fractional  turning  points,  such  equations 
admit  highly  irregular  points  at  which  the  coefficients  of  the  differential 
equation  can  be  almost  arbitrarily  multivalued.  Genuine  coalescence  of 
singular  points,  however,  is  not  considered.  A  local  representation  of 
solution  structure  is  established  which  generalizes  Frobenius'  method  of  power 
series.  Some  results  on  solution  symmetry  have  striking,  global  implications 
in  the  shortwave  limit . 
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SIGNIFICANCE  AND  EXPLANATION 


This  work  concerns  the  modulation  of  waves  or  oscillating  systems,  which 
pervade  all  the  science  and  engineering  disciplines  .  Modulation  occurs  when 
waves  travel  through  an  inhomogeneous  material  in  which  the  local  propagation 
velocity  differs  from  place  to  place,  but  the  differences  are  small  over  a 
distance  of  only  a  wavelength  —  a  very  common  case  in  the  sciences  and 
engineering.  The  resulting  change  to  the  waves  is  mostly  gradual,  but 
occasionally  drastic,  as  at  a  shadow-boundary,  where  oscillation  turns  into 
decay  and  quiescence  over  just  a  few  wavelengths  .  When  this  phenomenon  can  be 
analyzed  via  an  ordinary  differential  equation,  such  a  boundary  is  called  a 
transition  point. 

At  first,  only  the  simplest  transition  points  representing  the  most 
typical  shadow  boundaries  were  studied.  But  then  some  phenomena,  such  as  wave 
reflection  and  scattering  cross-sections,  came  to  be  traced  to  hidden 
transition  points  that  become  visible  only  when  real  distance  (or  time)  is 
embedded  in  its  complex  plane.  When  the  material  properties  vary  in  a  general 
manner,  (which  can  often  be  observed  only  incompletely)  the  hidden  transition 


points  can  have  arbitrarily  complex  structure  .  The  following  work  explores 


that  structure  in  detail  in  order  to  contribute  to  the  technical  basis  for  a 


reform  of  the  theory  that  will  make  it  simpler  and  more 
furnish  the  tools  for  more  efficient  scattering  calculations. 


IRREGULAR  POINTS  OP  MODULATION 


R.  E.  Meyer  and  J.  F.  Painter 

1 .  Introduction 

The  general  linear  second-order  equation 

e^d^w/dz^  +  q^w(z)  =0  (1) 

with  parameter  e  and  analytic  coefficient  function  q(z)  is  central  to  a 
vast  clast,  of  problems  in  the  sciences.  In  turn,  it  has  long  been  appreciated 
that  an  understanding  of  the  singular  points  and  turning  points  of  (1)  is 
central  to  the  treatment  of  those  problems,  regardless  of  whether  those  points 
occur  at  real  z.  Despite  many  studies,  no  complete  theory  of  such  points  has 
been  achieved,  perhaps,  because  (1)  encompasses  too  many  disparate  phenomena 
for  a  useful  theory  convering  them  all. 

The  present  study  focusses  on  those  forms  of  (1)  which  can  describe  the 
physical  modulation  of  waves  or  oscillators.  To  attempt  only  one  step  at  a 
time,  moreover,  it  excludes  questions  involving  genuine  coalescence  of 
singular  points  and  also  singular  points  of  (1)  artificially  introduced  as 
representations  of  radiation  conditions .  This  leaves  a  large  class  for  study, 
none  the  less,  because  the  coefficient  functions  q(z)  of  (1)  in  the  sciences 
must  be  defined,  if  not  by  speculation,  then  by  measurement,  in  which  case 
they  can  be  known  only  imperfectly.  The  characterization  of  q(z)  cannot 
therefore  be  very  specific  and  arbitrarily  irregular  points  must  be  admitted 
on  physical  grounds,  particularly  when  they  do  not  occur  at  the  real  values 
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of  z  of  direct  physical  meaning.  Certainly/  multivalued  functions  q(z) 
must  be  the  norm,  rather  than  the  exception.  All  the  same,  modulation  implies 
a  certain  structure  for  q(z)  (Section  2). 

The  main  objective  of  the  following  is  to  extend  to  the  general 
modulation  equation  the  insight  familiar  from  the  theory  of  regular  points 
that  the  local  structure  of  the  singularity  can  be  expressed  by  a  fundamental 
solution  pair  (wm,  wg)  of  (1)  associated  with  different  'exponents'  or 
'indices'.  Thus,  wm/ws  *  0  as  the  singular  point  is  approached,  so  that 
wm(z)  has  a  milder  singularity  than  wfl(z).  The  structure  of  irregular 
points  can  be  described  similarly  by  distinct  singularities,  one  stronger 
(Section  4)  and  the  other,  milder  (Section  3),  which  are  associated  with  a 
fundamental  solution  pair  constructed  below. 

In  the  general  case  of  modulation,  the  coefficient  function  q(z)  is  so 
poorly  specified  that  the  characterization  of  singularity  structure  cannot  be 
very  specific  either.  For  more  concrete  results,  the  general  class  of  such 
equations  is  therefore  slightly  restricted  (Section  2)  by  the  hypothesis  that 
a  limit  of 

q  2(dq/dz)  /*  q(t)dt 

be  identifiable  as  the  singular  point  z  =  0  is  approached.  This  still 
admits  a  much  larger  class  than  treated  before  [1]  and  makes  possible  a 
representation  extending  the  familiar  one  for  regular  points  in  terms  of  power 
series  commonly  associated  with  the  name  of  Frobenius  [2] .  The  representation 
established  below  is  in  terms  of  quasi-power-series  with  coefficients  that  are 
themselves  multivalued  functions.  The  results  reported  concentrate  on  certain 
aspects  of  those  functions  and  on  solution  symmetries  which  are  important  for 
the  applications  to  scattering  theory  that  motivate  the  present  investigation. 


A  clearer  indication  of  the  structure  thus  revealed  will  be  aided  by  a 
comment  on  the  role  of  the  parameter  e  in  ( 1 ) ,  which  is  proportional  to 
Planck's  h  in  the  quantum-mechnical  interpretation  of  that  equation.  The 
results  are  local,  in  the  first  place,  and  as  such,  concern  the  "parameter¬ 
less  case"  in  which  e  is  merely  a  fixed  complex  number  (say,  unity,  without 
loss  of  generality).  Since  all  the  proofs  are  on  that  basis,  e  could  in 
fact  be  regarded  as  excess  baggage.  Quite  the  opposite  is  true,  however, 
because  |c|  will  reveal  itself  as  a  homotopy  parameter  linking  irregular 
points  to  regular  ones  and  thereby,  as  a  signpost  to  a  broad  avenue  of 
approach  to  a  general  theory  of  irregular  points.  For  modulation  equations, 
moreover,  e  can  also  serve  as  an  illuminating  label  for  structural 
distinctions  even  in  the  parameter-less  case.  They  reflect  a  distinction  in 

the  scientific  meaning  of  the  two  terms  of  (1);  the  first  —  or  more 
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precisely,  w  d  w/d(z/e)  —  represents  the  oscillatory  mechanism  of  Newton's 

2 

law;  by  contrast,  q  (z)  represents  the  material  potential.  This  is  a  little 
too  naive,  more  rational  variables  x  and  ex  must  be  introduced  (Section  2) 
in  the  place  of  z/e  and  z,  respectively,  but  those  do  turn  out  to  play 
fundamentally  different  roles  of  oscillation  and  modulation  variables  even  in 
the  parameter-less  case.  A  third  reason  for  treasuring  c  is  that  the 
structural  results  have  a  striking  significance  for  the  shortwave  case 
I e |  <<  1. 

The  distinction  between  the  modulation  role  of  z  and  oscillation  role 
of  z/e  in  (1)  is  an  echo  of  the  distinction  between  'gravitational'  and 
'inertial'  mass  discussed  by  philosophers  of  physics  in  the  late  nineteenth 
century.  It  is  therefore  relevant  to  report  that  our  comment  on  the  benefits 
of  e  is  purely  mathematical  hindsight.  The  present  investigation  set  out 
merely  to  establish  semi -abstract  existence  constructions  adequate  to  support. 
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a  minimal  connection  theory  (3]  needed  for  scattering.  That  construction 
revealed  a  striking  'two-scale*  structure  in  terms  of  functions  of  x  and 
functions  of  £x.  The  functions  of  x  are  related  to  powers;  those  of  ex, 
by  contrast,  have  a  structure  peculiarly  different  from  powers,  they  are 
multivalued  'mild'  functions  (Section  2)  of  a  very  general  type  including 
logarithms  as  familiar  examples.  A  two-scale  structure  is  thereby  really 
implied  only  for  the  asymptotic  representation  (Section  5),  which  is 
remarkable  in  providing  information  global  in  x,  even  though  only  local  in 
ex.  Even  in  the  parameter-less  case,  however,  where  our  results  are  entirely 
local,  the  label  e  helps  to  recognize  the  distinction  between  the  role  of 
the  mild  functions  and  the  power  functions  in  the  solution  structure. 

For  regular  points,  certain  exponent  differences  lead  to  the  "Frobenius 
exceptions"  where  the  stronger  solution  contains  logarithmic  multiples  of  the 
milder  one.  Analogous  features  arise  for  irregular  points  (Section  4),  where 
they  are  less  pronounced  because  the  structure  is  already  pervaded  by  the 
'mild',  multivalued  functions.  However,  they  do  spoil  the  symmetry  of  the 
stronger  solution  (Section  5)  for  a  non-generic  subset  of  equations. 

One  would  naturally  wish  to  extend  to  general  irregular  points  Langer’s 
[1]  triumph  of  uniform  approximation  of  solutions.  It  will  indeed  emerge 
strikingly  below  how  general  modulation  equations  are  a  direct  extension  of 
Bessel's  equation.  Beyond  his  class  of  fractional  turning  points,  however, 
Bessel  functions  cannot  serve  as  uniform  approximands,  because  the  branch 
points  are  more  complicated,  and  it  is  not  readily  apparent  yet  how  far  the 
notion  of  uniform  approximation  could  be  extended  to  more  general  classes  in  a 
practically  useful  way.  To  some  degree,  at  least,  the  'local'  representations 
established  below  do  provide  such  an  extension  because  they  are  global  enough 
to  answer  asymptotic  questions  (3}. 


Good  error  bounds  for  the  approximations  proved  below  are  also  very 
desirable,  but  it  appears  premature  to  consider  them  in  detail  in  this  first 
foray.  There  are  many  indications  that  the  formulation  and  results  developed 
below  are  not  definitive,  but  merely  document  an  open  road  not  noticed 
before.  As  it  is  explored  further,  better  estimates  will  be  found.  Even  the 
present  ones,  however,  are  based  on  the  method  of  Volterra  equations  normally 
employed  [4]  to  obtain  error  bounds,  and  to  promote  access  to  them,  the  symbol 
6  will  be  used  to  denote  the  elements  from  which  such  bounds  wuld  have  to  be 
assembled. 

2.  Modulation  Equations 

Equation  (1)  is  only  one  of  a  family  of  normal  forms  of  the  general, 
linear,  second-order  differential  equation  and  constructive  statements  are 
awkward  and  cumbersome  in  such  an  indefinite  frame.  By  contrast,  the 
Lionville-Green  or  WKB  or  Langer  variable  x  such  that 

dx/dz  =  iq/e 

has  long  been  recognized  as  the  natural  one  for  the  description  of  waves  or 
oscillators;  it  measures  length  or  time  in  units  of  (2ir  times)  the  local 
wavelength  or  period,  and  physical  specifications,  e.q.,  radiation  conditions 
for  scattering,  relate  directly  to  x  or  ex.  Any  simple  and  practical 
theory  should  emerge  in  terms  of  these  variables,  and  continued  reference  to 
the  independent  variable  z  can  only  obscure  matters  and  will  be  abandoned 
after  this  Section. 

For  y(x)  =  w(z),  (1)  transforms  into 

y"  +  2fy'  =  y,  2f(x)  =  -ieq  2dq/dz  ,  (2) 

which  shows  modulation  to  be  controlled  by  the  function  f(x),  rather  than 
by  q(z)  directly.  It  explains  also  why  roots  of  q(z)  (’’turning  points”) 
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which  are  ordinary  points  of  (1)  play  the  same  role  in  the  theory  [4]  as 
singular  points  of  q (z)j  all  are  singular  points  of  the  modulation  function 
f(x)  and  hence,  are  the  singular  points  of  (2)  or  of  any  other  rational 
formulation  of  (1)  as  an  equation  describing  modulation  of  waves. 

Since  genuine  coalescence  of  singular  points  [5]  is  excluded  from 
consideration  here,  each  can  be  studied  separately.  This  does  not  exclude 
functions  q(z;£)  in  (1)  with  pairs  or  clusters  of  singular  points 
approaching  each  other  as  | e |  decreases,  provided  only  that  they  do  not  do 
so  too  rapidly  for  a  rescaling  [6]  to  exist  into  a  formulation  (1)  in  which 
they  remain  bounded  apart  independently  of  £.  In  order  not  to  overload  the 
presentation  here,  the  rescaling  [6]  is  presumed  a  priori,  and  the  remaining 
dependence  of  q  on  £  (which  is  then  immaterial  for  what  follows)  is 
ignored.  Attention  may  then  be  directed  to  a  single,  singular  point  with  an 
e-independent  neighborhood  N  free  of  others. 

To  define  the  subclass  of  equations  (1)  that  can  describe  wave 
modulation,  two  requirements  now  suffice.  First,  the  natural  variable  x 
must  be  defined,  for  otherwise,  not  even  the  concepts  of  wavelength  or  period 
could  exist  for  (1).  Secondly,  if  q(z)  be  non-integrable  at  a  singular 
point,  then  that  point  corresponds  to  no  £x  e  <C  and  hence  represents  not  a 
genuine  singularity  of  modulation,  but  a  device  for  re-interpreting  a 
radiation  condition  as  a  singular  point  in  the  z-plane.  This  will  be  excluded 
here  to  concentrate  on  the  class  of  genuine  modulation  equations.  For  it,  the 
singular  point  z  must  correspond  to  a  definite  point  £x,  and  without  loss 
of  generality,  both  may  be  identified  with  the  origin,  so  that 

Ex  -  i  q(t)dt  (3) 

exists  on  a  neighborhood  N'  C  N  of  z  *  0,  even  if  not  as  a  single-valued 
function,  and  this  is  the  main  premise  of  the  present  study. 
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To  make  it  effective  still  requires  transformation  to  the  framework  of 
the  natural  variables,  and  a  preliminary  remark  on  notation  may  be  helpful. 
Observe  that  (2)  shows  e  'f  to  be  a  function  of  z,  not  of  z/e,  and 
therefore  by  (3),  of  ex,  rather  than  of  x.  The  primes  in  (2),  of  course, 
denote  d/dx,  and  thus  y  must  be  anticipated  to  be  a  function  both  of  x 
and  of  the  parametric  variable 

Ex  =  5  . 

Explicitly  notation  to  that  effect  is  cumbersome  and  technically  redundant, 
because  the  proofs  below  are  for  fixed  e  ^  0,  but  perception  of  the 
structure  of  the  theory  will  be  made  easier  by  the  convention  that  a  notation 
such  as  q(x)  denoter  a  function  that  may  depend  also  on  5  while  a  notation 
such  as  ip(G)  denotes  a  function  (such  as  E  ^f)  only  of  the  modulation 
variable  5  =■  Ex. 

Now,  the  5-image  of  N'  will  contain  a  disc  about  5=0.  Since  the 
analysis  will  be  local  in  5,  no  generality  is  lost  by  subjecting  its 
radius  E  also  to  a  bound  E(Y)  specified  later  for  the  convenience  of  the 
estimates.  This  disc,  with  a  cut,  will  be  the  5-domain  A  of  the  analysis; 
the  corresponding  x-domain  D  is  the  cut  disc  of  radius  E/|e|.  An  intrinsic 

V, 

statement  of  our  premise  (3)  is  therefore  that  a  branch  r(x)  of  q *=  must  be 
definable  as  an  analytic  function  on  D  so  that 

idz/d5  =  r  2 

is  integrable  to  5*0. 

It  may  be  conjectured  that  this  physical  hypothesis  i3  sufficient  for  a 
theory  of  irregular  points  of  modulation.  To  obtain  more  concrete  results, 
however,  we  add  a  second,  lesser  hypothesis  that  the  function  xf(x)  = 
(x/r)dr/dx  which,  like  E  ^f,  depends  only  on  Ex,  has  a  limit, 
xf(,  ♦  Y  e  r  as  5  ♦  0,  uniformly  in  A  . 
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Accordingly 


xYp( ex) 


(4) 


and 

(£/P)dP/d£  =  xf(x)  -  Y  -  ♦  (£)  (5) 

is  analytic  on  A  and 

♦  (£)  ♦  0  as  £  ♦  0,  uniformly  in  A  .  (6) 

In  view  of  the  restrictions  to  be  placed  on  the  radius  E  of  A  for  the 
estimates,  no  further  loss  of  generality  arises  from  assuming  !♦(£)!  also  to 
be  bounded  on  A. 

To  interpret  the  second  hypothesis,  note  that,  apart  from  a  constant 
factor, 

P(£)  -  exp  /5  T“Vt)dT  (7) 

whence  it  follows  from  (6)  that  P(£)  is  analytic  on  A  and  varies  near 

£  »  0  more  slowly  than  any  nonzero  power  of  £, 

W  >  0,  ItV1!  ♦  0  as  £  ♦  0  .  (8) 

The  postulated  limit  Y  is  therefore  the  "nearest  power"  of  x  in  r(x)  * 

1/_  J/_ 

q  2  (and  q  2  is  recalled  to  be  the  familiar  amplitude  function  of  the  WKB- 
approximations  for  (1))*  The  integrability  premise  is  now  equivalent  to 

Re  Y  «  j  (9) 

with  an  added  restriction  on  P  and  $  in  case  Se  Y  ■  In  turn,  the  term 
"mild  function"  would  seem  apposite  for  any  function  ♦(£)  sharing  the 
defining  property  of  P(£)  that  it  is  analytic  on  A  and  (£/<|>)d<|'/d£  ♦  0  as 
£  ♦  0,  uniformly  in  A,  whence  also  (8)  then  follows  for 

For  further  interpretation,  note  that  the  singular  point  is  regular  for 
e  ■  0  because  then  $  =  0  and  the  modulation  function  f(x)  in  (2)  has  a 
simple  pole.  The  irregularity  function  ♦(ex)  defined  by  the  second 
hypothesis,  however  vaguely,  identifies  a  diffeomorphism  linking  the  irregular 
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points  to  regular  ones.  A  mathematical  role  of  homotopy  parameter  is  thereby 
revealed  explicitly  for  |e|. 

The  class  of  singular  points  of  Schroedinger  equations  admitted  by  the 
two  hypotheses  includes  very  irregular  ones,  due  to  the  unrestricted 
multivaluedness  of  $(£x)  and  hence,  also  of  r(x)  and  q(z),  in  addition 
to  all  the  turning  points  of  modulation  covered  in  the  literature.  It  extends 
even  the  class  of  [7]  (where  only  an  asymptotic  approximation  for  |x|  ® 

and  |e|  0  is  established)  by  abandonment  of  any  restriction  on  how  slowly 

$(?)  +  0  with  £.  For  Langer's  [1)  class  of  fractional  turning  points, 

2Y/( 2Y-1 ) 

z  q(z)  is  analytic  and  nonzero  at  z  -  0, 


<MC)  =  I  Yn  i 

n=  1 


;  cn( 1-2Y) 


and  the  solutions  of  (1)  and  (2)  are  approximable  in  terms  of  Bessel 
functions. 


3.  The  Milder  Solution 

The  regular  singular  points  of  (2)  have  Frobenius  exponents  or  indices 

[4,  p.  150]  0  and  1  -  2Y  and  the  auxiliary  function  z(x)  which  (3) 

associates  with  each  modulation  equation  (2)  will  turn  out  to  generalize  the 
1-2Y 

Frobenius  power  x  to  irregular  points.  Indeed,  by  (3)  to  (6)  and 

L'Hopital's  rule,  as  x  0  in  D, 

-2 

x  dz  _  -i£x  d(r  x)/dx 

2  dx  2  __  dz/dx 

rz  x+n  (10 

=*  lim  (1  -  2xf )  =  1  -  2Y 
x+0 


so  that 


and  £(S)  is  mild: 


1  -2  y 

z ( x)  =  x  C (  £x ) 


U/OdC/dC  -  «(C)  ♦  0  as  C  ♦  0 


uniformly  in  A. 

The  idea  behind  most  of  the  estimates  to  follow  is  that  such  mild 
functions  can  be  controlled  at  the  expense  of  arbitrarily  small  powers.  For 
instance,  if 


lub  |<Mu£)  +  u(uC)l  •  , 

ue(o,i) 

6 

P(u€)t(u£)u  V[P(CK(C)1  -  e^u,?) 

then 


( 12) 


lej  <  for  u  e  [0,1]  and  ^(C)  ♦  0  as  i  *  0  (13) 

because  <o  and  $  ♦  0  with  £  and 

log  e^u.5)  *  -/^  [^1  +  $(*£)  +  w(rO] 
has  non-positive  real  part.  A  restriction  |C|  <  with  an  E^  >  0  can 
therefore  assure  ^(C)  <  1  and  thereby,  6^  +  Re  Y  <  3/2,  by  (9). 

Theorem  1.  For  |€|  <  E(Y)  so  restricted  that 

M  =  (3/2)  -  Y  -  «,(£>  M4) 

has  positive  real  part  (and  at  worst,  |£|  <  E^),  (2)  possesses  a  solution 

y  (x)  *  z(x)y(x)  =  x1  2YC(tx)  [  a  (ex)(x/2)2n  (15) 


on 


D  with  z,  C  and  6  defined  by  (3),  (11)  and  (12)  and  or  given 


recursively  by  aQ  =  1  and 

fl  2n-1,  rl 


^(C)  *  4  /J  u2n  ’du  (e^  1,u5)  3  2an_^  (  ^uC)  ^2M  3  2ndl  . 
Proof.  By  (2)  and  (3),  y(x)/z(x)  -  y(x)  obeys 

y"  +  2 ( f  +  z'/z)y'  "  y  ■  0  , 


(16) 


which  will  be  satisfied  by  a  differentiable  solution  of  the  Volterra  equation 

$Mx)  -  f*  rliziiili.]2  $<v)dv 

y  J0  lr(x)z(x)J  y v vi av 

(17) 

y(x)  -  1  +  y’(v)dv  . 
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Let  m  =  Re  M  and 


^n+l 

dx 


2 

rx  r(v)z(  v). 
J0  r ( x) z ( x) 


a  (v)dv 
n 


a„  =  1,  a  (0)  »  0  for  n  >  1 
0  n 

then  by  (4),  (11)  and  (12), 


(18) 


-1 

x 


n+1 


ll 


1-Y-«1  2 

[e,(u,C)u  ]  a  (ux)du 

i  n 


(19) 


and  the  restriction  on  E(Y)  assures  by  (13)  that 

I  x  ^ajl  <  (2m)  \  a^(x)  =  -j  x2a^(£),  I  <»1 1  <  1/m  , 

and  if  a  (x)  ■  <*  (£)(x/2)2n  with  |a  (C)|  <  k  ,  that 
n  n  n  n 

| (x/2) ~2n_1a^+1 |  <  kn/(m+n),  an+1(x)  *  an+1 (C) (x/2) 2n+2 

|o„+1(5)l  <  VB+n)"1(n+1>’1  =  kn+1  * 


Inductively,  therefore 

a  (x)  =  a  (£) (x/2) 2n,  |a  (£)|  <  k  =  r (m) / [nl Hm+n) ]  (20) 

n  n  n  n 

and  these  bounds  show  the  series  £  a*(x)  of  functions  analytic  on  D  to 

n 

converge  uniformly  on  compact  subsets  of  D  to  a  function  y'(x)  analytic 
on  D.  The  series  £  a^  converges  similarly  to  y(x),  and  summation  of  (18) 
confirms  that  y(x)  satisfies  (17).  The  formula  for  oM£)  follows  from 
(18)  and  (19). 

Corollaries.  Since  z(0)  =  0,  by  (3),  also 

ym(0)  =  0  .  (21) 

For  e  =  0  and  Re  Y  <  tj#  in  which  case  the  singular  point  is  regular,  S 

and  a  are  constant  and  (15)  is  precisely  the  Frobenius  series  [4,  p.  149) 
n 

for  the  solution  of  (2)  with  the  property  (21).  (For  Re  Y  1  p  the 
integrability  premise  does  not  admit  t  *  0.)  For  e  /  0,  C(£)  and  o^(£) 
are  generally  multivalued  on  full  neighborhoods  of  the  irregular  point. 
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L'Hopital's  rule  shows  (Appendix  I)  none  the  less  that 

nia  (0)  -  -  Y)/r<n  +  |  -  Y>  . 

The  sense  in  which  y(x)  tends  to  an  even  function  of  x  as  £  ♦  0  will  be 
discussed  in  Section  5.  Summation  of  the  bounds  (20)  for  (15)  yields  growth 
bounds 

|y  (x)/z(x)|  -  |y(x)|  <  r(m) | x/2 | ’"“i  .(!x|) 

n  m-i 

|dy/dx|  <  r(a)|«/2|1~ai  (|x|)f  m-l-^-ReY-  «,(£) 

n  z  i 

for  the  milder  solution  ym(x)  in  terms  of  modified  Bessel  functions  [4, 
p.  60]  .  They  emphasize  even  further  how  the  irregular  points  of  modulation 
concern  generalization  of  Bessel  functions . 

4 .  The  Stronger  Solution 

Theorem  2.  If  -  Re  Y  is  not  a  positive  integer  and  E  is  suitably 
restricted,  then  (2)  has  a  solution 

m 

y  (X)  -  l  8  (ex)(x/2)2p 
8  0  p 

analytic  on  D,  with  6=1  and  6  (?)  defined  recursively  by  (23),  (29) 

Op 

and  Lemmas  1,2  below. 

Remarks.  Since  ya(0)  ■  1,  ys(x)  represents  a  stronger  solution.  So 

does  y8(x)  +  aym(x)  for  any  a  e  C,  but  the  clearest  representation  of  the 

branch  structure  of  the  singular  point  will  be  by  a  fundamental  system 

(y8,ym)  in  which  yg  is  free  of  such  additive  traces  of  the  branch  point 

of  ym.  Comparison  with  Theorem  1  shows  Theorem  2  to  complete  such  a  system 

because  6  (0)  exists  and  is  nonzero  (Appendix  II)  for  all  p  >  0.  In  fact, 
P 

l  8p(0)(x/2)2p  is  the  power  series  of  (28)  below  [4,  p.  60],  but  for 
t  ft  0,  the  6  (£)  are  generally  multivalued  on  a  full  neighborhood  of  the 


irregular  point.  Summation  of  the  bounds  of  Lemma  2  below  again  yields  growth 

bounds  for  y_  and  y'  in  terms  of  modified  Bessel  functions. 
s  s 

Proof.  A  twice  differentiable  solution  yg(x)  of 
r2y1 (x)  =  /*  r2y(v)dv  +  const. 

y(x)  =  1  +  /*  y1 (v) dv 

will  satisfy  (2).  For  Re  Y  <  -  4,  however,  the  normalization  to  y_(0)  =  1 

2  ® 

is  seen  by  (5)  to  require  a  regularization  of  the  first  integral,  and 
different  regularizations  will  add  different  multiples  of  the  milder  solution 
to  the  stronger.  A  procedure  avoiding  it  is  to  write  this  Volterra  equation 
as 

r2y' (x)  =  /*  r2y(v)dv  +  Cx  , 

cx  =  K  r2(y-yN)dv  (22) 


Y(x)  =  1  +  f*y’ (v)dv 

with  a  fixed  regularization  parameter  X  ^  0  and 


for 


N 


yN(x)  =  l  bp(x),  r  b^+1  =  /*  r  bp(v)dv 


y.-i  3  b-i  3  °'  b0  3 1f 


b  (x)  =  /*  b’  (v)  dv 

p+1  '0  p+1 


(23) 


This  implies 


0  <  p  <  N  <  -  -  -  Re  T  <  N+1 


'**»;>•  -  r\  - 


(24) 


(25) 


To  describe  the  structure  of  yN  for  N  >  1,  require  0  <  |ex|  <  E  to 

assure  definition  of  <)>(ex)  for  |x|  <  |X|  and  note  that 

lub  |$(S)|  =  <5  (  |  ex  | )  ♦  0  as  |  ex  |  ♦  0  , 

1C !< lex! 


(26) 


by  (6) ,  so  that  E  < 
outlines  a  proof  of 


with  an  Ej  >  0  will  assure  fi  <  1 .  Appendix  II 


Lemma  1.  For  |£|  <  E^  and  for  1  <  p  <  N, 

b  (x)  =■  8  { C )  ( x/2 ) 2p 
P  P 

and  |0  |  <  k'  ■  (k*  )p/pl  with  k*  (Y)  independent  of  x  and  X.  If 
P  P 

~  -  Re  Y  is  not  a  positive  integer,  the  result  remains  true  also  for  p  = 

N+1  for  sufficiently  small  E(Y)  >0. 

This  shows  yN(x)  to  be  a  ’polynomial'  in  x  of  degree  N  with 

coefficients  that  are  bounded,  but  generally  multivalued,  functions  of 

K  m  ex.  For  integer  -  ~  *  Re  Y,  the  bound  on  fails*  indeed 

(Appendix  II),  8  (0)  exists  for  0  <  p  <  N+1  when  -  -  -  Re  Y  is  not 
P  2 

integer,  but  only  for  0  <  p  <  N  when  it  is.  For  e  ■*  0,  in  fact,  y^(x) 
is  the  sum  of  the  first  N  terms  of  the  power  series  of  the  solution 

•1  -y 

lim  y  (x)  *  r(l  +  y)(x/2)2  I  ,  (x)  (28) 

e-0  8  2  y.l 

of  (2)  for  <M£)  5  0. 

To  attack  (22),  consider  now  a  sequence 


y  (x)  *  y„(x)  +  t  b  (x),  n  =  N+1,  N+2 , • • • 
n  N  M , ,  p 

p=N+1 

-  n  s{8»  V»*- 


(29) 


bp+1(0)  "  ° 


with  t^+1  given  by  Lemma  1  and  N  defined  by  (24),  so  that 
-  —  <  N+1  +  Re  Y  and  these  integrals  converge  for  p  >  N+1,  by  (4)  to  (6). 
Appendix  II  also  gives  a  proof  of 

Lemma  2.  If  -  |  -  Re  Y  >  0  and  not  integer  and  if  |C|  <  E(Y)  with 


sufficiently  small  E(Y)  >  0,  then 
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b  <x)  *  8  {£) (x/2) ' 
P  P 


pl|6  I  <  plk’  -  ( k* )  T(N+1  -  a)/r(p-s) 

P  P 

for  p  >  N+2,  with 

s  *  j  -  Re  Y  +  52(5)  >  0  , 

(30) 

6  (C)  «  lub  |$(u5)|  • 

ue(0,1) 

For  -  j  <  Re  Y  <  the  construction  simplifies  because  no 
regularization  is  needed'  Then  N  =  -1,  y^  =  0,  b^  =1  and  the  same  proof 
leads  to  Lemma  2  for  p  >  1 ,  except  that  s  <  0  and 

plk'  *  r(-a)/T(p-8) 

P 

The  sequence  (29)  of  functions  analytic  on  D  is  therefore  majorized  by 
the  partial  sums  of  a  power  series  convergent  for  all  x  and  converges 
uniformly  on  compact  subsets  of  D  to  a  function  ys(x)  analytic  on  D, 
which  summation  of  (29)  shows  to  satisfy  (22). 

For  the  sake  of  completeness,  information  on  the  nature  of  the  stronger 
solution  for  non-positive  integer  values  of  Re  Y  -  -  may  be  desired.  For 
e  =  0,  where  the  singular  point  is  regular,  these  are  the  exceptional 
exponent  differences  [4,  p.  150]  for  which  the  stronger  solution  has  a 
logrithmic  branch  point.  For  irregular  points  the  case  Y  m  •j  is  somewhat 
special  because  the  integrability  premise  does  not  then  admit  e  =  0,  but 
assures  ym(0)  ■  0»  by  contrast,  the  construction  of  the  stronger  solution 
makes  ya(0)  “  I  also  for  Re  Y  ■  j.  For  Re  Y  <  -  -j,  a  different 
representation  of  the  stronger  solution  can  help  to  elucidate  its  nature. 

Theorem  3.  If 


e[yjx>] 


2  /£  lr(v)]2bN(v)ym(v)dv  , 
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then  yN(x)  +  g(x)ym(x)  solves  (2).  For  non-integer  —  -  Re  Y  >  1,  the 
additional  specification 

P2x2Y-1[g<x)y  (x>  -  »m,1(0)(x/2)2N+2]  +0  as  x  ♦  0 
on 

makes  yN  +  gym  identical  with  the  stronger  solution  of  Theorem  2 . 

Proof.  For  any  twice  differentiable  function  g(x),  g(x)ym(x)  solves 
lr2(gym) *1 '  -  r2®.  ♦ 

because  ym(x)  satisfies  (2).  The  definition  of  dg/dx  -  g' (x)  by  the 

2  2  2 

theorem  and  (3)  makes  (r  y  g*  )*  *  r  ymt>N  and  by  (25),  therefore,  yN  +  gyffl 
solves  (2)  and  must  be  a  linear  combination  of  ym  and  ys .  But,  g(x)  can 
be  examined  with  the  help  of  Theorem  1  and  Lemma  1  (Appendix  II)s 
Lemma  3 .  For  non-integer  -  Re  Y  >  1 , 


-2N-2  ,  .  ,  .  .  .-4J-1 a  . 

x  g(x)ym(x)  ♦  4  &N+1(0)  as  K  ♦  0  , 

with  N  defined  by  (24)  and  8  (0),  by  (A7)  . 

N  +  i 

Since  N  >  0  and  yN(0)  *  bp  *  t,  it  follows  from  (21)  that  yN  +  gym 

represents  a  stronger  solution  y  of  (2)  normalized  to  y(0)  *  1.  From  (24), 

moreover,  2N+2  <  1  -  2Re  Y  <  2N+4,  while  Theorem  1  shows  ym(x)  *  0(x*  2Y?) , 

so  by  Lemma  3,  yN  +  gym  represents  the  same  solution  as  Theorem  2  up  to 
2)1+2 

terms  o(x  ).  The  choice  of  additive  constant  in  g(x)  specified  in 
Theorem  3  extends  this  agreement  sufficiently,  by  (10)  and  Theorem  1,  to 
preclude  a  difference  of  a  nonzero  multiple  of  ym(x)  between  the  stronger 
solutions  of  the  two  theorems . 

While  Lemma  3  holds  only  on  strips  of  it,  (31)  is  defined  by  Theorem  1 
and  Lemma  1  in  the  whole  left  half-plane  of  Re  Y  -  ^  and  there  depends 
analytically  on  Y,  by  (2)  and  (5),  when  x,  e  and  $(£)  are  fixed  so  that 
ym(x)  1*  0  •  In  the  strips,  the  proof  of  Lemma  3  (Appendix  II)  shows 


(x/g) dg/dx  ♦  211  +  2Y  +  1  -  2v  as  x  ♦  0  , 

and  as  Re  v  ♦  0  so  that  the  righthand  strip  boundary  is  approached,  g(x) 


is  seen  to  approach  a  function  varying  less  than  any  nonzero  power  of  x 
(although  it  approaches  a  strictly  mild  function  only  as  v  ♦  0, 

Y  -  j  ♦  integer).  Thus  g(x)  extends  the  logarithms  of  Frobenius  [4,  p.150] 
to  more  general,  irregular  points;  in  the  special  case  where 
P(£)  *  P(0)  ft  0,  as  for  a  fractional  turning  point  and  in  particular,  for  a 
regular  point,  equation  (All)  below  shows  g(x)  to  have  a  purely  logarithmic 
branch  point  when  Y  -  ~  is  a  negative  integer. 

5.  Symmetry  Bounds 

The  coefficient  of  (10)  depends  continuously  on  £,  at  any  fixed  x  ft  0, 
by  (3)  to  (6),  and  as  E  +  0,  the  equation  approaches  a  form  of  Bessel's  with 
solution 

lim  y(x)  =  T(|  -  Y)(x/2)Y-1/2l  (x) 

E=0  V2  -Y 

normalized  to  y(0)  =  1.  At  the  same  time,  the  restriction  on  |x| 

disappears  that  had  been  needed  to  respect  the  merely  local  definition  of 

$(C)»  80  that  D  expands  to  the  whole,  cut  plane  of  x.  As  £  ♦  0, 

therefore,  y(x)  =  y^z  tends  to  an  even  analytic  function  of  x  on  any 

compact  C  C.  The  same  is  true  of  y  (x),  by  (28),  as  long  as  —  -  Re  Y  is 

s  2 

not  an  integer.  It  is  of  considerable  interest  that  these  remarks  can  be 
sharpened  to  detailed  bounds  on  the  oddness  of  the  solutions  for  nonzero 
C  e  A,  which  can  serve  as  foundation  of  a  theory  of  asymptotic  connection 
[3]: 

Theorem  4.  For  x  and  xe  **  in  D  and  E(Y>  restricted  as  for 
Theorem  1 , 

ly(x)  -  y(  xe  1Ti)  |  <6  (  )  £ 1 )  T(m)  |  x/2 1  2  mi  (|x|) 

m  m 

and  5  ♦  0  as  |£l  *  0. 
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i 


Here  »  “  (3/2)  -  Re  Y  -  ^(C),  as  in  Theorem  1,  with  5^  defined  by 
(12),  and  Im  denotes  again  the  modified  Bessel  function  [4,  p.  60] >  a  proof 
is  given  in  Appendix  I.  In  case  the  relation  of  the  order  m  of  the  Bessel 
function  to  the  variable  5  be  considered  more  confusing  than  useful,  it  is 
of  course  permissible  to  proceed  from  (12)  with  6^  replaced  throughout  the 
proofs  by 

«#(Y)  -  luh  MO 

to  obtain  Theorem  1  with  M  ■  ( 3/2 j  -•  \  and  Theorem  4,  with 

m  *  3/2  -  Re  Y  - 

For  fixed  C  6  A,  the  oddness  of  y^/z  »  $(x)  can  therefore  grow  at 
most  exponentially  with  |x|. 

The  worst  part  of  our  story  is  a  proof  (Appendix  II)  of  a  similar  bound 
on  the  first  N  terms  of  the  stronger  solution  for  Re  Y  *  -3/2  (for 
Re  Y  >  -3/2,  yN(x)  is  strictly  even  because  y_1  s  0,  y  S  i): 

Lemma  4.  For  x  and  xe  **  in  D  and  sufficiently  restricted  E(Y), 

lyN(x)  -yN(xe"1,i)|  <  «*(KI)  Z  (k'(Y)|x2/4|]p/(p-1)I 

and  6'  ♦  0  as  |€|  ♦  0. 
s 

Finally,  the  bound  is  extended  in  Appendix  II  to  the  stronger  solutions 
Theorem  5.  For  non-integer  j  -  Re  Y(  x  and  xe  **  in  D  and 
sufficiently  restricted  E( Y) , 

ly.(x)  -y„(xe”Wi)|  <  c(Y)«fl( |5|)|x/2|2+8I  (|x|> 

Bo  "  *8 

and  *  0  as  Kl  ♦  0. 

Here  s  is  again  defined  by  (30),  but  the  same  proof  also  yields  the 

Theorem  a  fortiori  with 

8  »  -  ^  -  Re  Y  +  lub  } $(5)1  . 

2  5eA 

For  integer  j  •  Re  Y  >  0,  however.  Theorem  3  and  Lemma  3  show  Theorem  5 


to  fail 
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Appendix  I 


X(u,C)  -  i  {li<  |C!®i6)  -  u(|uC|ei6)}d6 
tend  to  zero  with  i,  so  that  also 

i  (|C|)  -  lub  lub  |1  -  exp  2X(u,t)|  ♦  0  as  I C I  ♦  0  . 

n  |t|<|C|  ue[0,1] 

Since  aQ  2  1,  <M>  gives  |x-1[a*(x)  +  a'(jx)]  I  <  j  fim/®  with 

m  -  Re  M  *  (3/2)  -  R»  Y  *  ^  ss  in  Theorem  1,  and  it  now  follows  recursively 

from  (Al )  and  (20)  that 

|(x/2)1"2nIa^(x)  +  a^(jx)]|  <  n25Bkn 

\(x/2)‘2n  [an(x)  -  a^(jx)]|  <  n  «B*n 
and  Theorem  4  follows  from  (20)  and  Theorem  1. 
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Appendix  II 


Proof  of  Lemma  1.  From  (23),  (4)  and  (5), 

b'(x)  =  r~2  /*  r2dv  =  x  /^/x  [e2(u,C)]2  u2Y+25du 

e2(u,C)  =  exp  [<J>(t)  -  6]  . 


(A2) 


Choose  arg  X  so  that  |arg  5  -  arg(ex)|  <  r  for  all  ?  e  A  (Figure  1)  and 
the  path  in  (A2),  so  that  uC  moves  on  an  arc  of  constant  modulus  from  ex 
to  5'  *  |ex|exp(i  arg  K)  (Figure  1)  and  thence,  radially  inward  to  £.  For 
u€  on  this  path  P, 


le2 |  <  exp(r6) 

with  <S  defined  by  (26)  .  For  fixed  5,  u  moves  also  along  a  path  partly  at 
constant  modulus  and  partly  radial,  and  since  N  >  1  confines  our  attention 
to  He  Y  <  -3/2,  the  restriction  |€|  <  E2  assures 

s1  =  Re  S1  >  0  for  S1  -  -Y  -  -  4  (A3) 


and  from  (A2) 


x  ’b'fx)  =  J*^X  (e2 (■£,£) ]2t  1  dt 

with  path  again  partly  radial  and  partly  around  (at  most  half  of)  the  unit 
circle,  and  it  follows  for  |x/X|  4  i  that 

|x  ^’(x)!  <  e2ir5  /g  It  1  dt|  <  ^  kjj  . 

By  (23),  therefore,  |4x  2b  (x) |  <  k' .  If  b  «(x)  =  8  (£)(x/2)2p  2, 

1  U  P  1  p—  I 

moreover ,  then  again  by  ( 23 ) , 


,./x  ,  2S  *1-2p 

-  2  f*/X  e,8  (S/t)t  dt 

1  Z  p“  I 

and  the  restriction  |5|  <  E2  assures  for  p  4  N  that 
2s1  +  1  -  2p  >  1  -  24  >  -1,  by  (24)  a. id  (A3).  Therefore 


,  i  ,  23  +  1-2p 

I  (2/x)  2p~'b’  (x)  |  <  2e2  lub  |3  J  /jt  d 

p  I C I < I cx |  P_1  0 


and  if  | 


3  |  <  k'  and 

P-1  p-1 


k’  -  2e 


1<n<N-1 


,  2s  -1-2n 

ll  It  ’  atl 


with  path  as  before,  then 

(2/x) 2pb  (x)  =  3  ( £ ) ,  13  |  <  k'k'  ,/p  =  k*  ,  (A4) 

P  P  P  P~'  P 

and  for  p  <  N,  the  Lemma  follows  by  induction. 

For  non-integer  j  -  Re  Y,  (24)  implies  N  +  Re  Y  <  -  so  by  (A3), 

2s^  -  1  -  2N  +  26  =  -2(N+1  +  Re  Y)  >  -1  and  for  fixed  Re  Y,  (26)  shows  an 

E(Y)  >  0  to  exist  which  makes  also  2s^  -  1  -  2N  >  -1.  The  integral  for 
2d-1 

(2/x)  *  b'(x)  then  exists  also  for  p  =  N+1  and  (A4)  remains  valid  with 

P 

k'  maximized  over  1  <  n  <  N. 

While  the  rough  bound  (A4)  suffices  for  present  purposes,  it 

underestimates  the  decrease  of  1 8  I  with  increasing  p.  A  more 

P 

representative  indication  of  it  is  obtained  by  computing  8p(0):  As  £  ♦  0, 

by  Lemma  1  and  L'Hopital's  rule, 

8  (£)  +  22p  1  p  1  lim  [x1  2pb"(x)] 

P  x+0  P 

if  that  limit  exists,  and  from  (23),  (5)  and  (6), 


1-2p.  ..  .  .  -2  1-2p  fx  2.  , 

x  b'(x)  =  (r  x  )  Jv  r  b  ,dv 
P  *  P“ ' 

♦  22“2p  lim  [3  ,(£)/(2p  -  1  +  2xr'/r) ] 

x>0  P_1 

-  2 1 -2p 8  .(0)/(p  +  Y  -  ~) 

p-1  2 

if  that  exists.  Since  8^  =  1,  it  follows  that 

8p(0)  =  p'Vp  +  Y  -J)'l8p-1<°)  * 

provided  Y  -  —  is  not  a  negative  integer,  and  then  by  induction 

pi  8p(0)  -  T(^  +  Y)/r<p  +  £  +  Y) 


for  all  p  >  0.  In  any  cape,  (24)  implies  N  +  j  +  Re  If  <  0  so  that  the 
induction  yields 

pi  8p(0)  -  <-1)Pr(j  -  y  -  p)/r(i  -  Y)  for  0  <  p  <  N  (A7) 

even  when  Y  -  ^  is  a  negative  integer,  but  8.(0)  is  then  seen  not  to 
exist. 

Proof  of  Lemma  2.  By  (6), 

$,<5)  <  5.  ■  lub  I ♦( C) I  +  0  as  E  +  0  (A8) 

Ced 

and  for  fixed  Re  Y  <  -  a  sufficient  restriction  of  E(Y)  assures  s  >  0. 
If  again 

Vi'*’  *  Vilt,lx/2|2P’2'  'Vi1  * 

for  some  p  >  N+2,  then  from  (29),  (4)  and  (7), 

2p-i  ,f  2  2(Y-«  +p-1) 

b*(x)  -  2(x/2)‘dp^1  Jl  (e,(u,C)rB  (uC)u  2  du  , 
p  03  p-1 

6 

e3(u,C)  -  u  2P(uC)/P(C)  -  exp  t *2( C>  +  ♦(tC)] 
and  by  (30) , 


|e3(u,C) I  <  1  for  u  e  (0,1]  .  (A9> 

For  p  >  N+2,  (24)  implies  p-1  >  -  j  -  Re  Y  and  therefore  also 

ReY-62+p-j-p-1  -  a  >  0  for  fixed  Y  and  sufficiently  small 
|€|  >0.  Hence, 

|(x/2)1'2pb*|  4  k*  ,/(p-1-s) 

p  p-1 

and  by  (29), 

I (2/x)  2pb  |  <  k*  /[p(p  -  1  -  s)]  -  k*  (A10) 

P  P1  p 

and  Lemma  2  follows  inductively  from  the  bound  for  b^+1  of  Lemma  1. 

Proof  of  Lemma  3.  By  Lemma  1  and  (3),  (4)  and  (15),  the  definition  (31) 
may  be  written,  with  N  +  j  +  y  *  v,  as 
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„  2“2N  r-2L_ )2  X-2N"2  r*  d*  8  $  v2-v1 

2  dx  ^  o  .  '  /n  ,.  v 


fjc  r  z  Q  A 

•'o  —  v v  av 


’  2'2N"'  I  Sr  V01  -S’10  38  « * 0  • 


(All) 


by  (10),  (17)  and  (A7).  Since  (24)  implies  -1  <  Re  v  <  0  and  since  P  is 
mild,  dg/dx  Is  seen  not  to  be  integrable  to  x  =  0  unless  Re  v  =  0  and 
P(C)  *  0  as  £  ♦  0. 

2  2  V 

For  Re  v  <  o,  L'Hopital's  rule  may  therefore  be  applied  to  p  x  /g 
and  yields,  by  (All),  (5)  and  (6), 

.2  2 V,  ...  2xp-1dg/dx  +  2V  ... 

/9  7ThT  -  ■  2v/cn 

P  x  dg/dx 

as  x  0,  and  by  (All),  (x/g)dg/dx  ♦  2v,  which  therefore  represents  the 
nearest  power  in  g(x).  It  then  follows  from  (4),  (10),  (11),  (A6)  and  (All) 


-2N-2  ,  t  ,  ,  .a  -2V 
x  g(x)y  (x)  =  Cyx  g 


ec 

N  _  -N-1. 

2iv( 1 -2Y)  4  n+ 1  0 


as  x  ♦  0. 


Proof  of  Lemma  4.  From  (23),  for  1  <  p  4  N  defined  by  (24),  and  again 
with  exp(-ui)  =  j, 

b.(v)  ♦  b-ov. .  rx  ©Virait  *  '(f  ♦  ^x»>f^,Vt,s,ds 


=  (v/2)2p_1(I1  +  I2  +  I3)  , 

i,(v,  -  (|)2p '  rx  O'V'"1  •  V',)t>)at  • 


fit  Ir(y)1  ll  *  ®  ■ 

Vv>  •  ^2!> '  /?  'fnTr,2V',l>at  • 

To  obtain  bounds  on  these  integrals,  note  from  (4)  and  (5)  that,  if  ev,  tuv, 


jev  and  jtuv  are  all  in  A, 
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r( juv)r(v)  rCv  $(T)-$( jt) 


.  ruuvmv)  _  rev 
9  r(uv)r(jv)  "  ^euv 


r1  $(  eyo)  -$(  euva) 


do  -  8(u,ev) 


say,  and  |$|  <  6(|cx|),  by  (26).  If  ev  and  euv  lie  on  the  path  P  of 
Figure  1  from  €  to  ex,  then  the  integral  may  be  evaluated  along  the 
appropriate  part  of  P  and  from  the  first  form  of  8,  the  circular  part  is 
seen  to  contribute  at  most  2x6  to  (8|.  From  the  second  form,  the 
contribution  of  the  radial  part  is  seen  to  have  the  same  bound,  and  together 


| 8 |  <  4x6  , 


6  ( | ex | >  -  lub  lub  |1  -  e |  *  0 

3  IC|<|ex|  P 


(A12) 


I  ex  I  <  0  . 


Moreover,  if  p(es)/p(jev)  ■  e  (es,ev)  is  considered  for  ev  e  &  and 

4 

s  »  | X | exp  io,  0  >  o  >  -xf  then  (Figure  1) 

/-|ev|exp  io  .  dX  /|ex|  .  io  dr 
109  ®4  “  *  jev  *(X)  T  +  '|ev!  >  7“ 

and  by  (30),  (26)  and  (6), 

|e4l  <  e  4  |x/v|°  , 

6  (K|)  ■  lub  6  (5)  ♦  0  as  |€|  *  0  . 

4  e  2 

arg  4 


(A13) 


From  Lemma  1,  b  (x)  -  b  _(jx)  =  6  . (?) (x/2) ,  and  if  it  be  now 

p-1  P~ i  P~ * 

supposed  that  1 8  _.,(€)  I  *  (p-1 )  X  _^k*  ^ ,  then  by  (4)  and  (A2), 

IXjl  *  2 1 / jj/v  [e2(u,ev)]2  6p_1  ( eUv)u2(Y+5+P  ^du| 

<  2.J,s(p-i)ip.1k^1  ;’/v  ,uJp'2s'-3au, 

with  Sj  defined  by  (A3),  so  that  again  Re  +  1  -  p  >  1  -  6  >  0,  by  (24) 
and  if  E( Y)  is  restricted  as  for  the  first  part  of  Leimaa  1.  As  in  the  proof 
of  that  Lemma,  therefore,  it  is  deduced  that 

I1,!  ‘  <P-,>1p-il<'  • 


I 


r 


Similarly,  by  (A12)  and  as  in  the  proof  of  Lemma  1,  if  k*  ■  lub  | S  |  as 

P  P 


there  defined, 

ll2l  -  21/’  e*(l  -  e2S(u'evJ }6  ( jeuv)u2’Y+fi+p”” 


'X/x 

<  k'J  k1  , 
3  p-1 


du  | 


With  path  at  constant  |t|  *  |X|, 


I3  -  2  /^V  j"2Yte4(et,ev)]2Bp_1(et)(t/v)2(Y+p"1)d(t/v) 

and  again  with  8^  «  Re  *  -Re  Y  -  ^  -  6, 

2ff5  2(s  +1-p) 

|13I  <  2ire  V_1  |v/X| 


Therefore 

3 

b  (x)  -  b  (jx)  ■  /*  I  I,<v)(^)  p~  dv  -  8  (£)(£)* 
P  P  0^12  p2 


with 


2ir6 


2(s  +1-p) 


IB  (5)1  -  21/’  [  I  (ux)u2p“1du| 

P  U  1  A 

<  l(p-i>V,  *  V  f  ki-i  *  vT  •" 'V<|x/X'"’’ 

2*6  2(s  +1-p) 

-  ^[(p-DXp.!  +  53  +  —  e  ^7  Ix/X|  ] 

by  (A4).  Since  Aq  ■  0  because  bQ(x)  =  1,  it  follows  recursively  for 
1  <  p  <  N  that 

b  (x)  -  b  (jx)  -  3  (C)(x/2)2p,  |8|  <  pA  k’ 

P  p  P  P  P  P 

2*6,  .  12(81+1-p>  ,  2 
*  *  r  2*  4  lx |  r.  .  p-1  ,x, 

\  <  S  +  k'  (31~  e  Ix|  [1  +  V  'x1 

n-2  x  ^  lx  2p-2 

♦*—111  ♦•••♦Jlfl  ] 

P  A  p  X 

where  it  should  be  recalled  from  (24)  and  Lemma  1  that 

0<1-6<s.+1-N<2-5  and  k'  -  lk'(Y)]P/p! 

1  P 


(A14) 


Finally, 


|yN(x>  -  YN(jx)l  -  1 1  8p(C)(|)  |  <  l  pApkpl'jl 


( A15) 
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To  prove  Lemma  4,  it  therefore  remains  to  establish  a  bound  6'  on  X 

3  p 


which  tends  to  zero  with  |€|.  For  sufficiently  small  |x/X|, 

4* 


V  V,EX,)  +  krTs^+TT  *x/x* 


2(3^1-?) 


( A16) 


by  (A13),  and  E(Y)  has  been  restricted  so  that  s1  +  1-  p>1-<5>0  for 

1  <  p  <  N,  with  6  defined  by  (26).  The  requirements  on  the  regularization 

parameter  |X|  used  so  far  are  that  it  be  independent  of  x  and  that 

| C|  <  |ex|  <  E(Y),  the  radius  of  the  cut  disc  A  on  which  the  estimates 

hold.  For  fixed  Y  and  E(Y),  envisage  now  a  sequence  of  subdiscs  about 

C  *  0  of  radius  |ex  |  <  E(Y)  and  tending  to  zero  as  n  ♦  ",  and  consider 

n 

any  point  sequence  {C^}  C  A  such  that  I C^ I  +  0  as  i  ♦  Given  n,  the 

n-th  subdisc  will  contain  all  C^  with  sufficiently  large  i  and  there  is  a 

subsequence  (n(i)}  wuch  that  n(i)  ♦  “  and  |x./X  |  +0  as  i  ♦  For 

i  n(  i) 

all  JCJ  <  JC^I#  regularization  can  be  achieved  with  the  choice  I  f°r 

| X | ,  and  the  bounds  (A12)  and  (A16)  then  imply  a  sequence  decreasing 

towards  zero  as  i  *  *  and  such  that  X  <  5  for  |C|  <  |C. I  and 

pi  i 

1  <  p  <  N.  A  monotone  positive  function  S^(|C|)  therefore  exists  such  that 

X  <  «'(|C|)  for  1  <  p  <  N  and  5’<|Cf)  ♦  0  as  |CI/E  0  .  (A17) 

p  s  s 

Proof  of  Theorem  5.  By  Lemma  1,  the  computation  of  bp(x)  for  p  <  N 
from  (23)  extends  to  p  *  N+1  when  ^  -  Re  Y  is  not  integer.  In  that  case, 
(24)  implies  -  ^  -  Re  Y  “  N+A  with  A(y)  >  0,  and  for  sufficiently  small 
| ex | ,  s^  -  N  =  A  -  5  >  o,  by  (26),  and  the  estimate  (A14)  then  extends  also 
to  p  ■  N+1.  For  p  >  N+2,  on  the  other  hand,  bp  is  found  from  (29j  and 
Lemma  2,  so  that  then 


b' (v)  +  b' ( jv) 
P  P 


,»»»>« 


again  with  j  •*  exp(-ti).  As  in  the  proof  of  Lemma  2, 
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r(t)/r(v)  =  (t/v) 


Y-« 


e3(t/v,ev)  , 

and  as  in  the  proof  of  Lemma  4, 

lo,  -  I)  do  -  S(t/v,ev)  , 

r(t)r( jv)  3  o 

and  if  v  and  jv  are  in  D  and  t/v  =  u  e  [0,1],  then  this  integral  may 
be  evaluated  at  constant  |o|  to  see  by  (30)  and  (6), 


|0(u,eV)  |  <  2it32(ev)  , 

( A18) 

lub  lub  |1  -  exp  20(u,ev)|  =  5H(|£|)  +  0  as  |£|  +  0 
|v|<|x|  ue[0,1]  s 

Now  suppose 

Vi(x)  -  Vi(jx)  =  Vi(^(x/2>2p"2'  'Vi'  <  (p“1)Vik' 

tr  r  v  p—  1 

which  is  known  for  p-1  =  N+1  from  the  extension  of  (A14),  then  since 


b^(v)  +  b^(jv)  ~  v  Jq  e^Cb^^tuv)  *  e29(u,ev) b^_^ ( juv) }u  2du 
and  by  Lemma  2,  for  p  >  N+2, 

1-2p  ,  ,  „  20  2(p-1+Y-«  ) 

(j)  (bMv)  +  bMjv)]  =  2  JQ  e3tSp_1  +  d  -  e  )  8p_1  }u  du  , 

where  again  p-^-S2+ReY=p-1  -s  >  0  for  fixed  Y  and  sufficiently 

small  Kl  >  0,  by  (24)  and  (30).  By  (A9),  (A10)  and  Lemma  2,  therefore 

|(v/2)1“2p[b'(v)  +  b*(jv)]|  <  pk'[(p-1)X  +  S"}  , 

P  P  P  P~  1  S 

b  (x)  -  b  (jx)  =  3  ( 5) ( x/2 ) 2p  , 

P  P  P 

IB  I  <  k'[(p-1)X  +  5-]  . 

P  P  p- 1  s 

It  now  follows  recursively  from  (A14),  (A17)  and  (A18)  that 

|3  |  <  5  (Kl)pk'  for  all  p  >  1,  with 
pa  p 

5  (|C|>  -  max[5'( |CI),«"( KIM  +0  as  Kl  >  0  . 

8  S  S 

From  Theorem  2,  therefore 

00 

|ya(x)  -  yg(jx)|  <  $g(KI)  l  p>M|x/2|2p 
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